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Abstract: We explore regions of parameter space in a simple exponential model 

A Q 

of the form V = Vq e m p that are allowed by observational constraints. We find 
that the level of fine tuning in these models is not different from more sophisticated 
models of dark energy. We study a transient regime where the parameter A has to 
be less than a/3 and the fixed point Qq = 1 has not been reached. All values of the 
parameter A that lead to this transient regime are permitted. We also point out that 
this model can accelerate the universe today even for A > \/2, leading to a halt of the 
present acceleration of the universe in the future thus avoiding the horizon problem. 
We conclude that this model can not be discarded by current observations. 
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1. Introduction 

The universe today is in a stage of accelerated expansion and we still do not know 
exactly what is the physical mechanism that drives this acceleration. Several alter- 
natives have been proposed, such as the old cosmological constant or the effect of a 
field whose potential energy dominates nowadays the energy of the universe [[1]. This 
last proposal is usually called quintessence f2|. 

One of the simplest models for quintessence is a scalar field rolling down an 
exponential potential |§. Exponential potentials are natural in theories with extra 
compact dimensions, such as Kaluza-Klein supergravity and superstring models [|]]. 

Recently, it was pointed out the possibility of having exponential quintessence 
with a temporary acceleration phase, avoiding the horizon problem that appears in 
the context of string theory associated with eternal acceleration ||. 

Exponential quintessence has been studied extensively in the literature |], [7| and 
usually discarded from limits on big-bang nucleosynthesis. This is the motivation 
for introducing more complex models with, e.g. two scalar fields H and double 
exponential potentials ||. However, the arguments against the simple exponential 
model are based on the assumption that a fixed point solution for the quintessence 
evolution is reached early on in the universe. Therefore, one has to analyse the case 
where the field have not yet reached this fixed point. This is the goal of this article. 
We find regions of parameter space for exponential quintessence that are allowed by 
observational constraints and examine the level of fine tuning required for realistic 
solutions. 
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2. Fixed point or fine tuning? 



We adopt the potential 



V(Q) = V e X ~v , (2.1) 



where Q is the scalar quintessence field, M p = (87rG)~ 1//2 is the reduced Planck mass 
and A and Vq are parameters of the potential. 

The equation of motion for the quintessence field is given by: 

Q + 3HQ = -^-, (2.2) 
dQ 

where dot denotes differentiation with respect to time. Friedmann's equation for a 
flat universe is: 

H 2 = ^ = ^ (pm + P R + \q 2 + V(Q)} , (2.3) 

where R is the scale factor and H is the Hubble parameter. For scalar fields, the 
energy density and pressure are respectively pQ = \Q 2 + V{Q) andpg = \Q 2 — V(Q). 
The radiation and matter densities pR and pu evolve as: 

Pm(r) + nHpM(R) = , (2.4) 

with n = 3(4) for non-relativistic (relativistic) matter. 
The acceleration R of the universe is given by 

where i = M,R and Q, u>m = 0, ojr = 1/3 and ojq = Pq/pq are the equation of 
state for matter, radiation and quintessence respectively. The universe is accelerating 
today if ujq < — |, assuming that quintessence is currently the dominant component 
of the energy density. The index denotes the present epoch. 

The dynamics of the quintessence field is obtained by solving the non-linear, 
coupled differential equations (|2.2j ) and ( |2.3| ) with appropriate initial conditions. 
This system of equations has stable fixed points (fp) for A < ^Jn with VLq = 1 and 
A > \fn with VLq = p- 0], where n depends on which background component is more 
important when the field reaches the fixed point: n = 3(4) for matter (radiation). 
The case A > y/n is the selftuning (or scaling) scalar field discussed by Ferreira and 
Joyce and by Copeland, Liddle and Wands 0. In this case, however, either the 
system reaches its non-trivial fixed point early on in the universe and the value of 
D,qo is too low or the quintessence field energy density contributes too much to the 
total energy density in the early universe and spoils big bang nucleosynthesis (BBN) 
predictions. Furthermore, in this selftuning case the universe is not accelerating 
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today, since when the fixed point is reached, ujq imitates the dominant background 
component equation of state (matter today) and consequently ojqq = ujm = @, [?]]. 
For these reasons, this solution can be discarded. 

The case of interest is given by A < y/3, in which the quintessence has not yet 
reached the fixed point regime today but will do it in the future. In other words, Qq 
is evolving from a small value (in order that quintessence does not spoil BBN) to its 
value at the fixed point, Qq = 1. In this case, this is not a selftuning solution, but 



a tracking [JTOj one, where ujq is also frozen with the different value 



4 P = ^ > ( 2 - 6 ) 

being able to accelerate the universe today. Other regimes cannot explain BBN 
and/or the present acceleration of the universe, when one is using simple exponential 
potential. 

In all quintessence models there arises the so-called cosmic coincidence problem: 
why is quintessence starting to dominate just now or, in other words, why is the value 
of quintessence energy density nearly equal to the matter energy density today? No 
quintessence model has solved this problem until now, since, in order to explain 
this coincidence, in general one adjusts the potential energy density to make the 
selftuning (or the tracking) field to freeze in a energy density that is of the order 
of the critical density today. In [ffH], for instance, two examples are considered: 
V(Q) = K^ + ^Q~ a and V(Q) = K A [exp(M p /Q) - 1], where M p is the Planck mass 
and K is a constant. For any given V(Q), there is a family of tracking solutions 
parameterized by K, whose value is fixed by the measured value of fi M0 . This means 
that the overall scale of the potential for all quintessence models has to be of the 
order of the present energy density. For this reason, K is a free parameter used to 
fit current observations. Whether the order of magnitude of the initial conditions 
and parameters can be obtained in reasonable physical models is a very important 
question that has not been fully answered in a satisfactory manner. 

Hence, because the cosmic coincidence problem has not been explained by any 
quintessence model, what one can do is to discuss the naturalness of the choice of 
initial conditions and parameters. Our goal in this article is exactly this: we study 
what is the size of the region in parameter space and how precisely adjusted these 
initial conditions and parameters must be in order to explain some observational 
constraints, that is, what amount of fine tuning is necessary in simple exponential 
when one requires that it has not yet reached its fixed point regime. 



3. Relevant equations 

We will follow closely the approach described in Cline ||. Instead of integrating in 
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time, we will use another variable u defined by: 

u = ln(l + z) = -In R/R Q , (3.1) 

where z is the redshift, and we take Rq = 1. We also define new dimensionless 
variables: ^ ^ 

^ " 73X g ' A " 3Mj(HZf Pi ' (3 ' 2) 

-I ~ if 

^ = 3M p 2 (# m ) 2 7 ' H = H™ ' ^ 

where if™ = 100/io km s _1 Mpc -1 is the currently measured Hubble parameter and 
therefore: 

3M*(H™) 2 = p™ = 1.88 x 1(T 29 h 2 g/cm 3 . (3.4) 

The dimensionless variables pi and V are given in terms of the critical density. In 
our numerical investigations, we adopt the following values, for ho = 0.7: 

Pmo = PJ ^ = 0.3 and p m = ^ = 4.3 x lO" 5 ^ 2 = 8.5 x 1(T 5 . (3.5) 

PcO PcO 

The variable H is determined by the numerical solution H of Friedmann's equation 
divided by ii m . 
In terms of the new variables, the potential is given by 

V(Q) = Vo , (3.6) 

and the differential equations read: 

H 2 Q" - (p R + \pu + SVj Q' + ^ = ; (3.7) 
H> = P» + P*+J ; (3.8) 

i - ±Q' 

where prime denotes differentiation with respect to u. This corrects a typographical 
error in Cline ||, where he writes V instead of 3V in equation ( |3.7| ). Notice that 
because the universe is flat, one is able to isolate H in Friedmann's equation. 

These are supplemented by the solutions for the evolution of matter and radiation 

Pm = Pmo e 3u , (3.9) 

and 

PR = Pro e 4u , (3.10) 

where prq an d Pmo are the radiation and matter energy densities today. The quin- 
tessence energy density is given by: 

p Q = Ih 2 Q' 2 + V . (3.11) 
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With these new variables, one can also write 



that always satisfy 



q r + n M + n Q = 1 . (3.13) 

From the equation above, one can see that the values of Qmo, Qro and Qqq (since 
the universe is flat) depends on the solution of Friedmann's equation. The values 
(|3.5|) used to solve the equations will correspond to the measured values of VLmo and 
n m on\y if H = H n . 

The equation of state for quintessence is given by: 

Q 2 -2V H 2 Q' 2 -2V 

UJQ = ~ = n • (3.14) 

4. Numerical Results 

4.1 Parameters 

As stated before, in all quintessence models there is an overall constant (Vq in our 
case) that is determined by the fact that the major contribution to the energy of the 
field today must come from the potential term, and that the energy density of quin- 
tessence is approximately equal the present measured critical density (more precisely, 
the critical energy density minus the matter plus radiation energy densities). This 
work does not aim to discuss the naturalness of such a choice of Vq. We simply follow 
the same approach used in all quintessence models to "solve" the cosmic coincidence 
problem. 

From equation (|3.14| ) one sees that the pressure-to-density ratio (equation of 



state) has a current negative value only if the major contribution for the energy 
density of the quintessence field comes from the potential term. The observational 
fact that Qq ~ implies that 

V(0) = V Q e -^° « « (l(T 3 eV) 4 , (4.1) 

where V(0) denotes the present value of the potential energy density. The value of 
Vq depends on the value of Qq. As long as the initial condition Q' i is small (and 
we will see that this is the case from the equipartition of energy) and the field has 
not yet reached the fixed point, one can see that Qq ~ Qi, as figure [I] illustrates. 
Therefore, the values of V and Q, are related by equation ( |4.1|) . 

In this sense, the only free parameter in this model is A. As mentioned above, 
we will be interested in the range < A < a/3, since this region is the only one able 
to explain all observational constraints. 
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Figure 1: Value of the field Q as function of u. For the initial conditions showed, the field 
has approximately the same features. The field remains almost constant until it reaches 
the fixed point regime, when it starts to roll-down the potential with a constant derivative 
in u. The initial conditions used along the text are showed by solid line. These curves were 
made with Vq = 1 and A = 1.3. 

4.2 Initial Conditions 

In order to solve the differential equations, two extra parameters are needed, namely 
the initial conditions Qi and Q\. We are interested in the case in which the field has 
not yet entered in the fixed point regime. In this case, the field has today almost 
the same value it had initially, and one is able to relate V^(0) and Qi, in such a way 
that Qi can be absorbed in the definition of Vq. For this reason, we will take, with 
no loss of generality, Qi = 0. We will discuss later the role of Qi in the solution. In 
fact, changing the value of Qi just corresponds to a rescaling of the problem. 

The freedom in the choice of Qi also comes from the fact that the potential 
energy does not contribute to the initial value of density parameter of the field, VLq^. 
This can be seen by noticing that ^(0) is of the order of the critical density today (see 
equation (|4.1| )) and that our initial conditions are taken at z ~ 10 13 . This implies 
that Vq is at least 39 orders smaller than the energy density scale (critical density) 
of that epoch, since p c oc H 2 oc R~ 3 (R~ 4 ) during matter (radiation) domination 
epoch. The initial energy of the field is then in the form of kinetic energy, or 




We have assumed a flat universe, which implies that \Q'\ < The equation above 
shows what would be needed in order to have the limit case \Q'\ = V2 : the dominant 
energy of the universe must be the kinetic energy of the field. 
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Natural initial conditions from equipartition of energy after inflation suggests 
that Q,Q t i ~ 10~ 3 [[U|. From equation (fO|) one has then Q\ « 0.05. 



Therefore, the most likely set of initial conditions are the values Qi = and 
Q' { = 0.05. Nevertheless, we will also study the effect of taking all the possible set 
of initial conditions, namely Qi > and < \Q\\ < y/2. 

4.3 Constraints 

We evaluate the equations and demand that the solutions must satisfy some observa- 
tional constraints. The first is given by nucleosynthesis. Nucleosynthesis predictions 



claims [II] that at 95% confidence level: 

fi Q (lMeV « z = 10 10 ) < 0.045 . (4.3) 

Another observational constraint is given by the quintessence density parameter. 
Observations from cosmic background radiation anisotropy indicate that the universe 



is flat. A set of complementary observations indicates that Hmo = 0.3 ± 0.1 flTJ 
Thus, 

n Q0 = 0.7 ±0.1 . (4.4) 

The uncertainty on Qmo implies that there is an uncertainty on the determination 
of Vq. This is the reason why one can study a region on parameter space (Vq, A) in 
spite of the fact that Vq is not a free parameter. 

The last observational constraint considered here is the present quintessence 



equation of state [12] 



-1 < ^qo < -0.6 . (4.5) 

The fact that the universe is accelerating today, as will be shown in figure f|, does 
not imply that it will accelerate forever. The present value of equation of state of 
quintessence in the cases studied in this work is temporary, since the equation of 
state is frozen only when the field reaches the fixed point regime, what will happen 
in near future. 

4.4 Results and Discussion 

We solve numerically the coupled differential equations (|3.7|,^.8[) using equations 



( P^B| , |3TP| , |3.10| ). The effect of all possible initial conditions will be discussed later. 
First, however, the particular choice (the most likely one) of initial conditions Qi = 
and Q\ = 0.05 will be studied. 

For this set of initial conditions, the region of parameter space 1 able to satisfy 
all observational constraints, contrary to common belief, is reasonable, varying from 
A = to A ~ 1.7 and from Vo ~ 0.5 to Vq ~ 3.0. In other words, all possible values 
of A in the tracking regime are able to satisfy the present observational constraints. 
The results are shown in figure 0. 

x It is useful to recall that we are interested in models with A < -\/3 ( see section [l]). 
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Figure 2: Region of parameter space that satisfies all observational constraints discussed 



in subsection |4.3| . There is a reasonable region of parameters of the exponential potential 
that can explain all observations. In fact, all values of A that produce the tracking solutions 
satisfy the constraints. The uncertainty on Vq comes from the uncertainty on f^A/o- 



Figures |3|, [| and || shows respectively how the energy densities, the equation of 
state of quintessence and density parameters varies with u. Initially, quintessence 
contributes to a small fraction of energy of the universe and decreases as i? -6 , dom- 
inated by the kinetic term (ujq = 1), faster than matter and radiation densities. 
When the potential term becomes important, there is a rapid change in the equation 
of state from ujq = 1 to ojq = — 1 and the quintessence density freezes until today, 
when it becomes dominant. Afterwards, the quintessence reaches the fixed point 
regime, which is characterised by the quintessence density parameter going from as 
small value at u « 2.5 and reaching the fixed point, Qq = 1, when u ~ —2.5. In the 
fixed point regime (tracking solution), the ratio between kinetic and potential energy 
densities becomes constant and consequently ojq is given by the equation (|2.(j|). In 

this regime, the energy density decreases as R~ 3( - 1+IjJ q ) = R~ x2 . 

Figure [| shows the behavior of equation of state for various parameters. Initially 
they have the same behavior, but they evolve differently after quintessence enters the 
fixed point regime, because different parameters correspond to different frozen ratios 
between kinetic and potential energy densities. A is the parameter that determines 
what value Uq will be frozen in: small values of A corresponds to small values of 
ujq and vice versa. In particular, for A > y/2 the universe stops to accelerate in the 
tracking solution, since uOq > — | in these cases. For small values (A < 0.15) the 
quintessence behaves like a cosmological constant in the tracking solution. 

The correlation between the free parameter A and the equation of state today is 
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Figure 3: Energy density as function of u for quintessence, matter and radiation for a 
typical solution. Initially quintessence density decreases faster than matter and radiation. 
When the potential term becomes important, pq freezes. When the fixed point regime 
(tracking solution) is reached, the field changes again its behavior and decreases as R~^ 2 , 
lower than matter. Plot made with parameters Vq = 1 and A = 1.3. 



shown in figure || Notice that there is a degeneracy for A > 0.6, namely, different 
values of A generate the same ujqq. This degeneracy is mainly due the fact that the 
present value of equation of state is not its value in the fixed point regime, since the 
tracking solution was not yet reached. In the fixed point regime the degeneracy does 
not exist, since equation ( |2.6| ) is satisfied. 

With better measurements of ujq , one could put severe constraints on the expo- 
nential potential model, specially if uq < —0.85. This is the region of low A and, as 
it was commented before, in this region quintessence behaves as a cosmological con- 
stant today. It is important to realize that the figure |^ remains the same for almost 
all sets of values of Vo, Qi and Q \ that are able to produce a tracking solution we are 
interested here 2 , since the value of equation of state depends mainly on A. In fact, 
in order to have this plot unchanged, it is enough that exists a considerable region 
of A in the parameter space (Vo,A). Later will be shown that the existence of this 
region only depends on Q' { . Only for values of Q\ > 1.25 this plot will be changed. 

Another important aspect is the dependence of these results on Qi. Figure |7| 
shows the region of parameter space that satisfies all observational constraints when 
one uses as initial conditions Qi = 1 and Q' \ = 0.05. When Qi ^ 0, the region of 
parameter space is deformed in V by a factor of e v ^ AQ,; , because in these models 
Vo can always be rescaled as -Do = Voe _v ^ AQl , since Q ~ Qi. Hence, in this sense, 

2 Qi > 0, < \Q'i\ < \/2 and Vq determined in the way discussed in subsection 4.1. 
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Figure 4: Equation of state as function of u for various parameters. Initially all have 
the same behavior. When quintessence reaches its fixed point, each solution is frozen in a 
specific value, depending mainly on A: small values of A correspond to smaller values of 
oj^q and vice versa. For a small A, quintessence behaves in the fixed point regime like a 
cosmological constant. For A > y/2, the universe does not accelerate in this regime. 



choosing another value of Qi corresponds to just a rescaling of the "old" region. 



As it was seen in equation ( |4.2|) , the initial condition Q' i corresponds to the initial 
value of density parameter of quintessence. The most likely value is Q \ = 0.05 
because of the energy equipartition after inflation. However, one can argue that 
choosing this value could also be a fine tuning. For this reason, it is important to 
explore what happens with the allowed region of the parameter space if Q' i has a 
different value, keeping in mind that < \Q'j\ < y/2. The results are shown in figure 



The region of parameter space becomes smaller in Vq when Q\ increases. This 
happens because the field has to satisfy the constraint ( f4.3p . When the density 
parameter of quintessence is large initially, in general the density parameter of quin- 
tessence contributes too much to the total energy density during the BBN epoch, and 
a smaller region of parameter space is able to satisfy the observational constraints. 
Nevertheless, a reasonable region of parameter space still exists for almost all values 
of Q'j. In fact, the region of parameter space only vanishes when > 0.75. In 
other words, for almost all possible values of Q 1 i in a flat universe there still is a 
significant region of the parameter space able to satisfy all observational constraints. 
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Figure 5: Density parameters of quintessence, matter and radiation. The fixed point 
regime is characterised by quintessence density parameter equal to unity (^q 1 = 1)- The 
fixed point solution was not yet reached today, which is a transition epoch. The parameters 
used were Vq = 1 and A = f .3. 
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Figure 6: Allowed region of (A, ojqq) space. Notice that ujqq < —0.6, according to con- 
straint (|4.5| ) . This plot is independent on values of Vo an d Qi and changes only for "high" 
values of Q\ {Q' i > 1.25). The fact that different values of A are able to give the same loqq 
comes from the fact that today the fixed point regime was not yet reached. 



5. Conclusion 



In this work we have studied the simplest quintessence model with an exponential 



- 11 - 



I 



Q\ = 0.05 
£l Q = 0.00125 



10 



20 
A 



30 



40 



vo 



Figure 7: Region of parameter space that satisfies all observational constraints, when 
one is using as initial conditions Qi = 1 and Q' i = 0.05. Note that using Qi / just 
corresponds to rescale the parameter space in Vq by a factor of e^ x ® 1 . 
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Figure 8: Regions of parameter space for various Q' ' i different from that of inflationary 
equipartition. Since the universe is flat, \Q\\ < \/2. For almost all possible values of \Q' i\ 
a significant region on parameter space still exists. It only vanishes for Qq^ > 0.75. 



potential. This potential is usually discarded because it cannot satisfy simultaneously 
all observational constraints. This is true in a regime where the field has already 
reached its fixed point regime. We have studied this model in a regime where the field 



has not yet reached its fixed point regime (tracking solution) today. We have shown 
that, contrary to common belief, this potential is able to satisfy all observational 
constraints for a reasonable region of parameter space (Vo, A). 

We have also shown that the resulting parameters and initial conditions are not 
unnatural. On the contrary, we showed that these parameters and initial conditions 
are not less natural than that used for other models of quintessence. For almost all 
possible values of Q\ there still is a significant allowed region of parameter space. V 
is determined by the measured value of Qmo, the same feature used in all quintessence 
models. Qi does not affect qualitatively the region of parameter space, since in the 
regime in which the field has not entered its tracking solution Qi and Vo are related, 
and consequently Qi just rescales the problem. The only free parameter of this model 
is A, which determines the behavior of the field. All values of A that produce the 
tracking solution satisfy the observational constraints. Depending on A, the universe 
may or may not accelerate forever. 

The allowed regions that we found are essentially due to present experimental 
uncertainties: the region of V is due the uncertainties on measured value of VL MQ and 
the region on A arises from the observational uncertainties on ujqq. If the uncertainties 
were reduced, these regions would also be reduced. In this case, one would be able 
to determine the values of the parameters of the simple exponential potential using 
observations. 

This potential cannot be discarded by any of the constraints discussed here, even 
if better measurements were made. The only way to discard this potential based on 
the constraints discussed in this work would be if constraint of BBN ( f4.3|) had been 
more stringent. For example, if one had showed that fig(lMeV) < 10~ 6 , probably 
the region of parameter space would become smaller, in such way that could not be 
possible to explain the measured value of u;qo, for example. Another way to discard 
this potential is to include another constraint as, for instance, the value of Qq during 
the epoch of formation of structure or in the last scattering surface, and verify that 
this model is not able to satisfy all constraints simultaneously. Therefore, we showed 
that at the moment there is no reason to discard the exponential potential or to 
consider it less natural than any other quintessence model. 
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